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APPLICATION OF THE METHOD OF CHARACTERISTICS TO SUPERSONIC ROTATIONAL FLOW 

By Antonio Fjskbi 


SUMMARY 

A system for calculating the physical properties oj supersonic 
rotational flow with axial symmetry and supersonic rotational 
flow in a two-dimensional field was determined by use oj the 
characteristics method. The system was applied to the study 
oj external and internal flow jar supersonic inlets with axial 
symmetry. For a circular conical inlet the shock that occurred 
at the lip oj the inlet became stronger as it approached the axis 
oj the inlet and became a normal shock at the axis. The region 
in which strong shock occurred increased with the increase oj 
the angle oj internal cone at the lip oj the inlet. For an inlet 
with a central body the method oj characteristics was applied to 
the design oj an internal-channel shape that, theoretically, results 
in very efficient recompression in the inlet; it was shown that 
ij an ejjuser is connected with the diffiuser a body oj revolution 
with very small shock-wave drag can be determined. 

INTRODUCTION 

The characteristics metliod for the determination of super- 
sonic phenomena was first used by Prandtl and Busemann 
for two-dimensional fiow (references 1 and 2). For flow 
with axial symmetry Frankl (reference 3) used the method 
of characteristics for determining tlie shape of a supersonic 
circular effuser with uniform exit velocity, and Ferrari 
(references 4 and 5) independently used the characteristics 
method for determining supersonic phenomena for every 
type of boundary condition. Subsequently Guderley (refer- 
ence 6) and Sauer (reference 7) transformed the system pro- 
posed by Frankl and Ferrari and obtained a different 
analytical solution of the problem. In all applications the 
hypothesis of potential flow was made; therefore the 
equation of potential flow was used. 

"When shock waves that are not plane (two-dimensional 
flow) or conical (flow with axial symmetry) occur in uniform 
flow, the variation of entropy across the shock is not constant 
and the flow behind the shock is no longer isentropic and 
becomes rotational. If the variation of entropy is small, 
the effect of rotation of the flow is not important for de- 
tennining the pressure distribution along a body and the 
theory of potential flow gives correct results. If the shock 
wave is strong and has large curvature, however, the effect 


of the rotation becomes important and the flow must be 
considered as rotational. 

The method of characteristics can be extended to apply to 
rotational flow if, in place of the potential function for the 
differential equation of motion, the stream function con- 
sidered by Crocco (reference 8) is used. With the character- 
istics method for rotational flow a more exact determination 
can be made of the shape of the shock wave and the distribu^ 
tion of velocity and pressure for phenomena in which the 
effect of rotation is important, as in the internal flow through 
supersonic inlets. The procedure of numerical calciilation is 
similar and not much more complicated than that used for 
the case of potential flow with axial symmetry. 
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SYxMBOLS 


pressure 

density 

entropy, mechanical units 

velocity 

!Mach nmnber 

limiting velocity corresponding to adiabatic expansion 
to zero pressure 
speed of sound 


ratio of velocity to 


limiting velocity 



ar-component of relative velocity 
y-component of relative velocity 
Cartesian coordinates 

Mach angle ^arc sin 

angle between velocity I’' and r-axis, radians 
angle between tangent to shock and direction of 
velocity of flow in front of shock 
deviation of direction of velocity across shock wave 
angle of polar coordinate in conical field 
potential fimction 

stream function for rotational flow (see equation (11)) 
ratio of specific heat at constant pressure and constant 
volume 
gas constant 
normal to streamline 
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H, 

fi = 
&= 

A= 

fif= 


sin /3 tan sin </> __ 

cos (vj+iS) 

sin |8 tan )3 sin y 
~ cos (<o— j8) 

L, K, and N defined by equation (14) 

I 

=Xa — Xb 
= Xc — Xa 

\Va 

"Wb 

=1+^ 


e = 


sin Pa 


f 

r= 

A 


cos (Pa-^-iPa) 
sin Pb 


cos (<pb — Pb) 
tan PA+h tan Pb 
Sb — Sa 1 

yii e+gj 


Subscripts: 

0 cbamber condition (zero-velocity adiabatic transfor- 
mation) 

A points of first family 

B points of second family 

0 quantities in the points calculated from A and B 
X derivative with respect to x 

y derivative with respect to y 

a ahead of shock 

b behind shock 

<pCn value corresponding to value of <p at point C» 

C„ at point C„ 

CHARACTERISTICS METHOD FOR SUPERSONIC POTENTIAL 
FLOW WITH AXIAL SYMMETRY 

The differential equation for potential flow with axial 
symmetry (reference 4) is 


\ a? ) 'by^ a* bx^'y ^ ' 

In supersonic flow some lines can be individuated (charac- 
teristic lines) that divide the flow into two regions for which 
the values of <t>, 4>x, and along the fine are different. For 
eveiy point of the flow two characteristic lines can be de- 
termined; every line is inclined at the Mach angle with 
respect to the direction of the velocity at the point, and 
therefore the characteristic lines can be divided into two 
families on the basis of the sign of the angle of the characteris- 
tic line with respect to the direction of the velocity. A 
family that is usually called the fii’st family is defined by 
tlie equation 

^=tan (^-h^) 


( 2 ) 


and the other family (second family) is defined by 
^=tan i<p-p) 


( 3 ) 


The variation of the quantities that define tlie velocity 
{({> and TO along a characteristic Une is given by the follow ng 
equations from reference 5: 

For the first family, 

dV dx 

^-d^tan;S-Z^=0 (4) 

and for the second family, 

tan jS— m ^=0 (5) 

y y 

where I and m are trigonometric expressions defined as 

j sin p sin </> tan p ' 

~ COS{P+,p) 

(G) 

sin g sm y tan P 
cos (,<p—p) 

If the direction of the velocity and tlie Alach number at 
two points near each other (points A and B in fig. 1) are 
known, the direction of the velocity and the Mach number 
at a point C given by the intersection of the two chametor- 
istic lines of different family starting from points A and B 
can be determined. Because the distances BC and AC are 
small, all the coefficients of equations (4) and (5) can bo 
considered constant and coincident to the corresponding 
values at points B and A, and tlie tangents to the character- 
istic lines at points A and B can be substituted for the dmr- 
acteristic lines from point A to point 0 and from point B 
to point 0. In this case the lines AC and BC arc straight 
lines. The line BC is inclined at an angle (p—P with respect 
to the K-axis and the line AC is inclined at an angle <p+P with 



Figube 1.— Oeometrical construction for determining point O b; tlte cbaractcrlstics method. 
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respect to the i-axis. Equation (4) can be applied for the 
line AO whore 

dx=Xc—Xjt_ 

y=VA 

V=V^ 

l=h 

and equation (5) can be applied for the line BC where 

dx=Xc—XB 

y=yn 

V=Vb 

/3=|3 b 

m=niB 

For practical use equations (4) and (5) are combined and 
transformed by means of the ratio TT of the velocity V to 
the limiting velocity This ratio is defined as 


,r) 


= 1 -1 r sin’ )3 

7—1 


( 7 ) 


and the following equations are obtained: 

/ Jir \ 

d<pi [t&n tan J=i— tan j3js 




Ia w. 


ya 

^^=tan Pa d<pA- 


{Xc—Xa)Ia 

Va 


Wc=WA + dWA 




( 8 ) 


( 9 ) 


( 10 ) 


With the method of characteristics (reference 4) it can 
be shown that, if a deviation of a streamline which wets the 
body occurs, the phenomena on the corner are regulated by 
the same laws that regulate the two-dimensional flow; 
therefore, the tangents of every characteristic line starting 
from the comer are known. If the initial flow conditions 
are known, the step-by-step calculation of aU the physical 
properties of the flow in the entire field is permitted, partic- 
ularly the calculation of the shape of the shock wave and 
the pressure distribution along any body of revolution 
with axial symmetrical flow in cases in which the hypothesis 
of potential flow is correct. 

CHARACTERISTICS METHOD FOR SUPERSONIC ROTATIONAL 
FLOW WITH AXIAL SYMMETRY 

Supersonic perfect flow is rotational when the flow is 
preceded by a shock wave and when the variation of entropy 
across the shock wave changes from point to point behind 
the shock. In this case the transformation of the fluid 
along every streamline is isentropic until another shock wave 
occurs in the fluid (reference 8). 


If a stream function i' is assumed to be defined by the 
following equations from reference 8: 


the equation 


4'r=yu{l—W^)'^^ 

1 

>Ar=-2/»(l-Tr ’)-^‘ } 


M) = 


jf(l — 


( 11 ) 


( 12 ) 


is a function of only the stream function f (reference 8); 
and therefore from the equations of state, continuity, energy, 
and steady motion the foUomng equation can be obtained: 


Equation (13) is a Monge-Amp^re equation, and if 


(13) 


H=1 

K=- 


Vl 

c* 

uv 


(14) 

two characteristic families with the following equations can 
be obtained: For the first family, 

I K^ 

dx~H V-ff* ^ 

and for the second family, 

dy K. jK^ 
dx H 

If n is the normal to the streamline, equation (11) 3delds 

2 


(15) 

(16) 

(17) 

(18) 


and 


curl VxV 1 j 1 cfs 
-=;^grad s=:^^ 


a‘ yj 

because s is constant along every streamline; therefore 

1 


1 1 ds 
2y Edn 


(19) 

( 20 '’ 

( 21 ) 


yW{l-W^)^ 
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and the folloiving expressions can be obtained: 


7—1 


\kvx=yUx(.l — W^) TF*)r 

Wx 


1_ 

-1 


y (22) 


When equations (21) and (22) are substituted in equations 
(15) and (17) and (16) and (18) and the Mach angle and the 
velocity are expressed in polar coordinates, equations (15) 
and (17) become . 


^=tan (/3+^) 


(23) 


(24) 


(first family) 

(second family) 
and equations (16) and (18). become 


S“tan 


dx j.dx ds sin®j3 

dW , . oj dx dxds sin® )3 „ 

^+tanffdy--m-g^;^-— ^--^-^°^0(secondfamily)(26) 

Equations (23) and (24) are identical to equations for 
potential flow (2) and (3), and equations (25) and (26) are 
similar to equations for potential flow (4) and (5), differing 
only by the terms that contain da. Equations (23) to (26) 
permit a step-by-step calculation of tlie entropy, intensity 
of velocity, and direction of the velocity if the initial and 
boundary conditions are known. If all physical properties 
are known for two points A and B and if the two points are 
close to each other, the tangents to the characteristic lines 
at the points A and B can be substituted for the characteristic 
lines with close approximation. In this way a point C can be 
determined as the intersection of the second characteristic 
line of point B with the first characteristic line of point A 
(fig. 1), because <p and |8 are known for the points A and B. 

For the characteristic lines of the first family, equation 
(25) gives the variation of tp and W from point A to point C, 
and all the coeflBcients are known and correspond to the 

cZs 

coefiBcients for point A ; only the term ^ is unknown. From 

equation (26) the variation of <p and W from point B to point 

da 

C can be determined, and all other terms, except -j-, are 


dn’ 

known and equal to the values for point B. The tei-m 


ds 


can be determined from the value of the entropy for points 
A and B. 

ds 

From figure 1, in equation (25) the term ^ can be writ- 
ten as 

d& _.. (ac— S a) cos {^a+<pa) 
dn (r< 7 — sin i5^ 


(27) 


and in equation (26) can be written as 

da X^B—Sg) cos (yg— fffl) 


dn 


{xc—Xb) sin |8 b 


(28) 


If the points A and B are close to each other and the varia- 
tion of entropy is gradual, equation (27) caij bo written as 


da 


8b— i 




da 


(Xc—Xj) 


sm 


cos {^a+<Pa) 


■(xc—Xb) 


sm /gfl 

cos {(Pb—^b) 


(29) 


and^ can be considered equal in equations (25) and (26), 

• ds T 

In this case ^ is known; therefore <p and IF can be calculated 

for point C. For practical calculations equations (25) and 
(26) should be transformed into two equations (equations (31) 
and (32)) each of which contains only one of the unknown 
terms dw and dtp. 

For simplicity, let 

' ^.=Xa—Xb 
^ a=Xg — Xa 

"='+1 

A(p = (pA — IpB 

Wa 


^^Wb 

sin Pa 

cos (/3^-f ¥>.i) 
sin Pb 

■' cos {ipB—M 
r=tan pB+h tan ^a 

A Sb — Sa 1 

yR e+gf 

Then the following equations can be obtained: 

rdvA=l—h—Aipi&nfiB+A{gJ sin* ^B+he sin* jS^) 

I io^A 


(30) 


Vb 


Va 


tan ^A dtp A — Aa sin* ^a + 
Sg=aA-{-AeyR 

<PC=^ <PA-\- dtpA 

Bc=BA+tnix 

7—1 


sin* j8c 


2 

yrisQ 


(bv 0 


1 

/ 

^0^1 A 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 


In equation (31) the terms and become very im- 

Vb Va 

portant as y approaches zero, near the axis of the body; 
therefore, in this region the distance between the points con- 
sidered must be reduced to obtain sufficient accuracy. 
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CHARACTERISTICS METHOD FOR TWO-DIMENSIONAL 
SUPERSONIC ROTATIONAL FLOW 

In two-dimensional flow an equation similar to equation 
(13) can be obtained if a special stream function defined by 
the following equation is assumed (reference 8): 

In this case the equation of motion (equation (13)) becomes 
(reference 8) 

-(i-Tr>)53(^-i)/»)=o (38) 

Equation (38), like equation (13), is a Monge-Amp^re 
equation and permits the determination of two equations that 
define the characteristic lines and two equations that give 
the variation of the velocity along the characteristic lines. 
The equation with transformations analogous to the case of 
three-dimensional flow can be written in the following form: 


^=tan i^+<p) 

(first family) 

(39) 

^=tan i<p—P) 

(second family) 

(40) 

c?lF , „ j , dx ds sin’ jS 

TT ^ ysm cos {p+<py 

= 0 (first family) 

(41) 


^ +tan d<p-^ ^ cos"^-ff) =Q (second family) (42) 


Equations for tp and dTT similar to equations (31) and (32) 
can be obtained from equations (39) to (42) by using equation 
d s 

(29) for the expression ^ as follows: 

T d<pA=l—h—ii<p tan j3s+.A(a/ sin* sin* /3^) (43) 


dW 

^j^=tan jS^ d<p—Ae sin* Pa 

(44) 

Sc = Sa -{-'AeyR 

(45) 

<Pc=<PA-\-d<pA 

(46) 

Wc=WA + dWA 

(47) 

sm Pc- 2 (ttv 0 

(48) 


DETERMINATION OF SHAPE OF SHOCK AND PRESSURE 
DISTRIBUTION ALONG A BODY 

The physical properties of supersonic flow past a body of 
revolution in axially symmetrical flow can be determined 
step by step by the use of equations (31) to (36). The S3re- 



Ficuee 2.— Practical aystem of calcolatiiiK tbe flo v field for a polntcd-nose body of rerolutlon. 


tern of calculation is similar to that for potential flow. If the^ 
body begins with a point (reference 5), a cone tangent to the 
body (%. 2) can be substituted for the front part of the body. 
If point A is the point at which the body can no longer be 
considered coincident with the cone, the velocity at point A 
is known from the cone calculations (references 9 and 10); 
therefore, the shape of the characteristic line AB of the first 
family can be designed, because tp and |S corresponding to 
different angles 6 are known from the cone calculations. 
At point A the body turns through an angle Atp and the flow 
undergoes a transformation that is determinable by the laws 
of two-dimensional flow. The velocity and direction of the 
flow after the deviation A<p, therefore, and the tangents to 
the new characteristic lines of the first family at point A can 
be designed. At a point B, near point A, the intensity and 
direction of the velocity are known; consequently, the tan- 
gent to the characteristic Hne of the second family can be 
designed at point B, and the point Ci can be located. At 
point Cl the physical properties can be calculated with the 
equations of potential fiow (equations (8) and (9)) because 
the shock in front of the body is conical. With the same 
system the point O* is determined from the point B„ on 
the shock wave. In order to determine the flow of the point 
C,+i on the shock wave, the equation across the shock and 
the equation for the characteristic lines of the first family 
must be used. If e is the angle between the tangent to the 
shock and the direction of velocity of flow in front of the 
shock, 5 the deviation of the direction of velocity across the 
shock wave, and the subscripts a and b denote the conditions 
ahead of and behind the shock, respectively, the shock equa- 
tions can be written in the following form (reference 11): 


tan e _ 2 F 1 . 7— 1 ~] 

tan (e— 5)~7-{-1 sin* (e— 5)”^ 2 J 

_1 /' y+l A/.* A. 

tan 5 \ 2 A/a* sin* e — 1 ^ e 


(49) 

(50) 
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Pi. 


Pa T + 


^ <U.‘ 


Pa_ 


lO 


sm'^ £- 


2t 


7— 


0 

0 


Pi 7+lVW sin“ e"^ 2 






(61) 

(62) 

(53) 

(54) 


Equations (51) to (53) show that the general condition for 
obtaining potential flow behind the shock when the flow in 
front of the shock is potential is that Ma sin e must be con- 
stant. Since the value of Ma in front of the shock is known 

from equations (49) to (54), the values of and 

can be determined as fimctions of 5 and therefore the equa- 
tion for the characteristic lines of the flrst family (equation 
(25)) becomes * . . 




sin® /3c, 


_dW dv> 
dS 11 c. 


-tan /3c„ Ic^ 


sm®/3cr ,d(sj-s„),"l „ 

where is the velocity behind the shock, s^^^ is the 

entropy, and — are coefficients that correspond to 

a deviation across the shock for which the velocity assumes 
the direction of <pca> H'c, is the velocity, and Sc, is the entropy 
at the point C*. When the value of rf^has been determined, 
the deviation across the shock 5 and the corresponding values 
of W and e can be determined at point 0,^.1. 

In order to determine the velocity on the body at a point 
Ai, equation (26) is used. At point Ai the value of <p is 
known because the direction of the flaw is tangent to the 

c? s 

body and, therefore, dtp is known; ^ is also known because 

in equation (27) the value of Sc corresponds to the value of 
s at the point Aj and is equal to the value at point A, which 
is known from the cone calculation. From the value of IF, 
the pressure on the body relative to the pressure p<! (pressure 
for zero; velocity from isentropic transformation from the 
conditions behind the shock) can be determined as 


^7=(1-W^*)^ (56) 

On the lip of the nose of an open-nose body a shock wave 
tliat is a two-dimensional shock occui’s; therefore, the tan- 
gent to the shock on the lip is known and the pressure and 
velocity behind the shock are also known (equations (49) 
to (54)). In figure 3 the line AB and the velocity behind 
the shock at point B are known. With equation (26) the 
velocity at point C can be calculated with the system just 
described, and from the point B the point D along the charac- 
teristic line of the first family can be calculated by equation 


(55). From point D the point E can bo calculated and a 
point F can be interpolated, whicli permits tlie determination 
of points G and H. Because tlie curvature of the shock 
near the lip is large, a point C that is very near point A must 
be selected so that correct interpolated values may bo ob- 
tained. The point E can be recalculated from point G, and 
the points I and M can be recalculated from point F in order 
to have a second approximation. 

For the case of two-dimensional flow the procedure is the 
same as. that for three-dimensional flow; equations (43) to 
(48) are used instead of the corresponding equations (cqup.- 
tions (31) to (36)). 

EXAMPLES OF APPLICATION OF CHARACTERISTICS 
METHOD FOR ROTATIONAL FLOW 

DETERMINATION OF SHOCK SHAPE AND PRESSURE DISTRIBUTION 

ALONG THE EXTERNAL SURFACE OF A SLENDER OPEN-NOSE BODY OF 

REVOLUTION 

Theoretical and experimental calculations were made to 
determine the shock shape and pressure distribution along 
the external surface of a slender open-nose body of revolu- 
tion. The body considered is the same body for which cal- 
culations of the external pressure distribution were made by 
Brown and Parker of the Langley Memorial Aeronautical 
Laboratory by use of the small-disturbance theory. The 
calculations were made for a free-stream Mach number of 
1.525, for which a schlicren photograpli taken during tests 
was available for comparison of the calculated and test 
results,.... 

In order to determine the importance of rotation of the 
flow, calculations for the front part of the body were also 
made with the potential-flow characteristic equations. The 
results of these calculations differed only slightly from those 
obtained, from the characteristics method for rotational flow 
because the curvature of tlie shock was very small. The 



FiQCEE 3.— Practical system ol calculating the flow field for a slender open^^osc body of 

revolution. 
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Figubb 4.— Practfcal system of calculating the eitemal flow around a slender open-nose body 
of revolution for .\f=l,626. 


calculations were begun witb tbe detennination of the two- 
dimensional shock on the lip of the nose. The practical 
system of the calculations is shown in figure 4, and in figure 
5 the calculated streamlines and shock-wave shape are com- 
pared with the diock-wave shape obtained from teat results. 
In figure 6 the pressure distribution calculated by the char- 
acteristics method is compared with the pressure distribution 
dete rmin ed by the small-disturbance theory. The small- 
disturbance theory undervalues the increase in pressure that 
occurs through the shock, but the differences in the results 
obtained by this method and those obtained by the charac- 
teristics method are small. 


o Experimental 



Figcss S. — Calculated streamlines and shock wave Tor a slender oi>en-nose bodf o{ revolution 
for 1.52S, showing experimental shock-wave shape. 


DETERMINATION OF SHOCK SHAPE, STBEAIHLINES, AND PRESSURE DIS- 
TRIBUTION ALONG THE INTERNAL SURFACE OF A SLENDER OPEN-NOSE 
BODY OF REVOLUTION 

Three slender open-nose bodies with different conical inlet 
angles are considered (figs. 7 to 12) and the supersonic part 
of the internal fiow is studied for a free-stream Mach number 
of 1.6. For this type of body the internal shock produced 
at the lip of the inlet has a very laige curvature and the effect 
of rotation is therefore very important. The calculations are 
extended to the region in which the Mach number is 1.0. 



Figubb 6 . — Pressnredlstributlonaloiigtbeexternalsurfftceofasleiidoropeii-iiosebody ofrevidutioii for Af"l,B25. (Data for amall-dlsturbaDca theory from work by Browa and Paiftar 

of the Langley Laboratory.) 
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Fioube 7.— Practical system of oaieulatlog the Internal supersonic-flaw Quantities for a 
slender open-nose body of revolntion for Af-1.6. 



Fiquee 9.~FractIcaI system of calculating the Intmial-supctsonic-flaw properties for a 
slender open-nose body of revolution for Af=1.6. 



Fiocbe S.— Calculated streamlines aind shock wave for Internal flow in a slender open-nose 
body of revolution for JW= 1.0. 



Fioube 10.— Calculated streamlines and shock wave tor Internal flow In a slender open-nose 
body of revolution for Af=1.6. 


APPLICATION OP THE METHOD OF CHARACTEEISITCS TO SUPEBSONIC ROTATIONAL FLOW 


119 



Fwobe 11.— Practical sjatem of calculating tbe tntemal superaonlc-flow propcrtfea Cor a 
slender open-nose bod? of revolution for Af-=1.6. 

The hypothesis is made that subsonic boundary conditions 
and stability considerations permit a subsonic flow at the 
end of the supersonic flow such as results from the calcula- 
tions. The results of the calculations show that at the 
axis of the inlet a normal shock always occurs and that the 
region in which a strong shock occurs (with subsonic velocity 
behind the shock) increases in extension with the increase 
in angle of the internal cone. When the cone angle ap- 
proaches zero and the shock is a Mach wave, a complete 
reflection occurs at the axis of the inlet and the extension of 
the strong shock is zero. 

For large internal cone angles (figs. 9 to 12) the shock is 
a simple shock that becomes normal in the central part of 
the body of revolution. After the shock the compression 
continues but the characteristic lines cannot form an enve- 
lope. For small internal cone angles compression occurs 
gradually and an envelope of Mach lines occurs. The shock 
therefore reflects from the central part and another shock is 
generated. The form of the reflected shock is shown in figiue 
8. The ratio of the diameter of the region in which a strong 
shock occurs to the diameter of the inlet as a function of the 
internal cone angle is shown for M=1.6 in figure 13. The 
results are interestmg for the practical design of supersonic 
inlets of slender shape because they show that for large cone 
angles the central part of the body of revolution, in which 
the compression is not very efficient, is large. 



FicrsE 12.— Calculated atieamlliies and shock ivave for Internal flow In a slender open-nose 
body ol revolution for .V=1.6. 

DETEBMINATION OF PHYSICAL PBOPBRTIES OF THE INTERNAL FLOW 
THROUGH AN INLET WITH A CENTRAL BODY 

For a Mach number of 1 .6 an analysis of the shape of an 
inlet with a central body was made to aid in obtaining high 
efficiency. Theoretically a supersonic diffuser with or with- 
out a central body and having no shock losses or shock drag 
can be obtained (reference 12); but for practical use it is 
convenient to accept small shock drag in order to avoid large 
friction drag. The inlet considered (fig. 14) has a 10° central 
cone. The deviation across the conical shock is 48'. The 
shock is reflected by a cylinder that forms the external part 
of the inlet. The reflected shock produces rotational flow 
behind the shock and the deviation across this shock on the 
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Fiauu 14. — Practical system of calculating the shape of the central body, of the streamlines, and of shock-vravo shape for a supersonic Inlet at Af— 1.6. 
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Fiodse is.— P ractical sjrstem ot cakolatlag the shape of the central body and of the streamlines of the tall of a body without pressure drag for Af~l.ft. 


cylinder is 1° and on the central body, 2° 43'. If the cen- 
tral body behind the shock has the same direction as the 
velocity, the shock vriU not be reflected and isentropic 
compression can be obtained behind the shock (flg. 14). 
The design of the central body, therefore, was deter- 
mined from the calculation of the corresponding streamline. 
The variation of the velocity along the external cylinder and 
the value of the exit velocity were fixed, and from this con- 
dition and the condition dependent on the shock, the ve- 
locity at every point was calculated. In order to avoid 
errors stream tubes were designed that permitted, on the 
basis of the ratio of the area in the region of uniform velocity 
to the area at the end of the stream tube, a check on the 
precision of the numerical 'calculations. The !NIach num- 
ber in the minimum section of the inlet was fixed at a value 
larger than 1.0 so that disturbances from the subsonic part 
of the flow would not cause instability. The value chosen 
was 1.09. 

If an effuser is connected with the diffuser, a body of 
revolution with very low shock drag can be obtained (fig. 
15). The only pressure losses are the losses across the two 
shocks, which are very small; but for practical applications 
the friction losses are larger than for the iatemal body alone. 
A balance of the pressure losses and friction losses must 
therefore be made in order to examine the possibility of 
practical use of this arrangement. 

CONCLUSIONS 

A system for calculating the physical properties of super- 
sonic rotational flow with axial symmetry and supersonic 
rotational flow in a two-dimensional field was determined 


by use of the characteristics method. Practical use of the 
system is based on a step-by-step procedure, which requires 
long numerical calculations; but the calculations for three- 
dimensional flow are of the same type as for potential flow 
and, therefore, can be used for the practical problems in 
which rotation is important. Some applications were made 
to determine the external and internal flow on bodies of 
revolution with axial symmetry, and the foUowTng con- 
clusions were indicated: 

1. The effect of rotation is not very important if the 
variation of entropy is small but is important in the study 
of internal flow, for which the variation of entropy is usually 
large. 

2. When the inlet is a circular conical channel, a shock is 
produced at the lip of the inlet that becomes stronger when 
the shock approaches the axis of the inlet and becomes a 
normal shock at the axis. The region in which a strong 
shock occurs (with subsonic velocity behind the shock) in- 
creases with the increase of the angle of internal cone. 

3. If an inlet with a central body is considered, the method 
of characteristics permits the design of an internal-channel 
shape that, theoretically, results in very efficient recom- 
pression in the inlet ; and, if an effuser is cormected with the 
diffuser, a body of revolution with very small shock-wave 
drag can be determined. 


Langley MEnoRLiL Aeronautical Labor.\tobt, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., April 29, 1946. 
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